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We study the statistical and dynamic properties of the systems characterized by an ultrametric 
space of states and translationary non-invariant symmetric transition matrices of the Parisi type 
subjected to "locally constant" randomization. Using the explicit expression for eigenvalues of such 
matrices, we compute the spectral density for the Gaussian distribution of matrix elements. We also 
' compute the averaged "survival probability" (SP) having sense of the probability to find a system 

in the initial state by time t. Using the similarity between the averaged SP for locally constant 
randomized Parisi matrices and the partition function of directed polymers on disordered trees, we 
' show that for times t > t cr (where t CI is some critical time) a "lacunary" structure of the ultrametric 

D , space occurs with the probability 1 — const/t. This means that the escape from some bounded areas 

^""^ . of the ultrametric space of states is locked and the kinetics is confined in these areas for infinitely 

long time. 

On 

The explosion of the interest to the statistical properties of ensembles of random matrices around 1950s has been 
Ch ' motivated, in first turn, by a series of physical problems in nuclear physics. Later the physical applications of the 
\ random matrix theory (RMT) were extended to mesoscopics and statistical physics. One can mention recent works 
[H, HJ, which deal with physical applications of the eigenvalue statistics of ensembles of random matrices. The typical 
problems of RMT concern the evaluation of the averaged spectral density, as well as distributions of " level spacings" 
under the supposition that all matrix elements are independent randomly distributed entries which take values in 
some specific ensemble Q. 
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Despite the standard RMT can be applied to a wide range of physical phenomena, RMT, in its basic form, does 
not cover a particular class of hierarchical complex systems important for glasses, real networks and proteins. The 
description of such complex systems deals with the concept of hierarchical organization of energy landscapes [1, [|| . A 
complex system is assumed to have a large number of metastable states corresponding to local minima in the potential 
energy landscape. With respect to the transition rates, the minima are suggested to be clustered in hierarchically 
nested basins of minima, i.e. larger basins consist of smaller basins, each of these consists of even smaller ones, etc. The 
basins of local energy minima are separated by a hierarchically arranged set of barriers: large basins are separated by 
£f) ' high barriers, and smaller basins within each larger one are separated by lower barriers. Since the transitions between 
the basins arc determined by the passages over the highest barriers separating them, the transitions between any two 
' local minima obey the "strong triangle inequality" for ultrametric distances. Hence, the hierarchy of basins possesses 
the ultrametric geometry which has the natural visualization by the uniform p-adic (i.e. p + 1-branching) Cayley 



I/"") ■ tree. 

qq ' It might be said that the systems possessing the ultrametric geometry are hardly compatible with the systems 
, possessing the Euclidean geometry. Recently, however, the interrelation between hyperbolic (i.e. tree-like) and 
k/ ' Euclidean geometries has become a topic of physics of disordered systems. It has been shown in Q that one can 
. £^ | construct random Gaussian translationinvariant potential landscape which reproduces all the essential features of the 
■ Parisi landscapes directly in any finite-dimensional Euclidian space. Also the question of isometric embedding of a 
uniform Cayley tree into the 3D Euclidean space has been discussed in the work [7| . The concept of ultrametricity has 
been implemented in a number of toy models describing transport phenomena and so-called "basin-to-basin" kinetics 
of disordered complex systems [1, 0, [l(|. These models deal with various aspects of the ultrametric diffusion - a 
certain type of stochastic motion in a nonrandom energy landscape with regular hierarchy encoded in the symmetric 
Parisi transition matrix, T. The model of ultrametric diffusion has been successfully applied to the dynamics of 
proteins fill ]. The eigenvalues, {Ay}, of the matrix T determine a hierarchy of relaxation times of the entire system 
and, hence, define the kinetics constrained to an ultrametric landscape. The values of Xj are well known [gl llll [l3j: 

r 

Xj = -ptfW - (i -p- 1 ) p^'t^" 1 (1) 

7'=7+l 
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were Ao = by definition, T^ 7 > is the element of the Parisi matrix T related to the hierarchical level 7 of the p-adic 
tree, and the summation runs up to the maximal level max[7] = T. Since T is a kinetic matrix, the sum of all matrix 
elements in each column of T is equal to zero. 

The structural randomness is inherent for any disordered system. Thus, the representation of the energy landscape 
by a regular hierarchy is rather ambiguous. In some cases, like in spin glass models, one supposes that the distribution 
of energy barriers is fully random. However, for ultrametric energy landscapes it is more reasonable to suppose that 
the barriers, belonging to the same hierarchical level, coincide only on the order of magnitude. So, it is naturally to 
think that the transition rates between basins of the same hierarchical level have random distribution around definite 
mean values prescribed by the ultrametric hierarchy. Besides that, in real experiments we deal usually not with 
individual molecule, but with an ensemble of molecules in the sample. Therefore, the computation of an observable is 
connected with the averaging of a dynamics in an individual ultrametric landscape over the ensemble of its realizations 
in the sample. 

To specify the model, consider two basins, belonging to the same hierarchal level. The dynamic trajectories leading 
from one basin to the other one, pass through the common "saddle point". The transition rates between the states in 
these basins depend only on the location of the " saddle point" , but not on the details of the behavior of the system 
in the basins. This scenario implies that all elements of the Parisi-type matrix are always equal inside a block, but 
they can fluctuate from one block to another one. Thus we arrive at the randomized "locally-constant Parisi (LCP) 
matrix" , T, it is depicted in Fig[TJt. Recall that T is the kinetic matrix, hence the sum of matrix elements in each 
column is zero. 




Each matrix element, Tij = q\ n , defines the probability to pass over the corresponding barrier. Hence, the 
dimcnsionless barrier energy (the "barrier height"), = -lng^ } {u { " ] > 0) can be introduced. As it has been 
pointed above, we suppose that 

uW = (u< n >) +£<**> (2) 

where (^u^ 1 ^ is the mean barrier height in a given hierarchy level, 7, and describes the fluctuations around this 

mean value. The mean values ^it^"'^ fix the structure of the nonrandom (i.e. regular) Parisi matrix. If the barrier 

heights grow (in average) linearly with 7, we set e - ^ ') = (qj 1 ^ = p~( Q+1 - )7 , where a is a linear landscape "slope" 
(a > 0) @,[ll|]. 

Thus the problem of stochastic motion in the randomized ultrametric landscapes has appeared. The first important 
question concerns the computation of the spectral density of the system described by the randomized Parisi-type 
symmetric matrix T shown in Fig[T^,. Besides that, we analyze some probabilistic properties of the averaged "survival 
probability" using the analogy with statistics of polymers on disordered trees. For any distribution of £y the spectral 
density of the matrix T can be, in principle, computed numerically. However since the locally-constant Parisi matrix 
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T defines the new class of random matrices not yet considered in the literature, it is worth understanding the analytic 

£y . In what follows we suppose that |£ 7 ™ 



structure of spectral density of T for the Gaussian distribution of £ 7 ™' ) . In what follows we suppose that |£ 7 ™^| <C 1 in 



([2|). That allow us to truncate the power series for g 7 in £ 7 at the linear term: 

gW«e-H n) )(l-4»)+0((4-)) 2 )) (3) 

The non-negativity of the matrix elements of T is ensured by the condition |£ 7 | <C 1. 

The randomized matrix T can be viewed as an ultrametrically organized system of barriers defined by a nonuniform 
p-adic tree - see FigJTJi (p = 2). Each matrix element of T represents a transition probability over the largest barrier 
separating any two different states located on the lowest level 7 = 0. For example, the states A and B are separated 
by the barrier in the branching point D, and the transition probability from A to B is q% , the states A and C are 
separated by the barrier in the branching point G and the transition probability from A to C is q% (see FigJIJi). 

The eigenvalues of LCP matrix have been derived for the first time in [l2j using the elements of the p-adic analysis 
[l3| . The construction of A 7 . n has a very transparent geometric interpretation. Define a pair of numbers (7, n), were 7 
is the hierarchy level (1 < 7 < T = max[7]) and n enumerates the blocks in the same hierarchy level (1 < n < p r_7 ), 
as it is shown in FigQJi. Then the eigenvalue A 7) « of LCP matrix for p-adic tree, reads: 

r 

A 7 .„ = -P 7 <Z 7 " ) -(1- P - 1 ) Yl P 7 '^ W 

7'— 7+1 

s 

In Eq.Q <7 7 are the elements of the matrix T shown in FigQJi, while the summation E (despite it resembles much 
(fTj)) contains some ambiguities. It is clear from FigQJi that any pair (7, n) fixes some vertex of the non-regular p-adic 
tree. The sum E in (|4|), contributing to A 7) n, runs from the point (7', n'), which is the upper nearest neighbor of the 
point (7, n) along a p-adic tree, towards the root point (r, 1). For example, let us compute the eigenvalue A 7= i jn= 3, 
p = 2 (see FigQJi). The first term in (|4]) is just the weighted contribution from the point E(7 = l,n = 3), while 
the sum E in this case is the weighted sum of two terms coming from the points F and G (see FigQJi). Hence, 



A lj3 = -2M 3; -(l-2- 1 )[2 2 gf +2*qP 

Let -P(£ 7 ™' ) ) be the distribution function of the fluctuations ^ 7 "' ) of matrix elements q\ (see ([3])). The spectral 
density, p(X) of the matrix T can be computed in the standard way: 

p(A) = i^(i(A-A 7 ,„)) (5) 

{7,"} 

where JV = p T is the total number of eigenvalues of the matrix T and (...) means averaging with the distribution 
function P(^ n) ). 

Suppose that the distribution of £ 7 is Gaussian and does not depend neither on the hierarchy level, 7, nor on the 
index n counting different blocks in the same hierarchy level: 

Speaking formally, since the fluctuations of £ 7 ™' ) in ([6]) arc not restricted, the distribution P(^ 7 ™' ) ) in Eq.© is inconsis- 
tent with the condition |£y | <C 1 introduced above. Recall that the last condition ensures the matrix elements g 7 ™^ 
in ((3|) to be non-negative. However in what follows we suppose (|(£ 7 ™') 2 ^ = o 2 <C 1, derive the spectral density and 

check our results for selfconsistency by comparing analytic and numeric results. With (j4|), ([5]) and (|6|) in hands we can 
easily compute the spectral density, p(A). Proceed first with an auxiliary computation of Q(A, 7) = (<5( A — A 7)n= i)) 

for n = 1 and arbitrary 7 (1 < 7 < T), where the averaging is carried out with f (Cn™ -1 "*)- The straightforward 
computations give the following result for Q(A, 7) at a <IC 1: 

f°° dr r°° dr f r 2 1 

0(A ' 7) = / x 2^ ^ < e ^(er 11 )^^^)} w / x 2^ cxp \ ix [A + ^ (r)] - T a2 ^( r )| ( ? ) 
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where 

r r 

v^T)=p-^ + (i- P - 1 ) J2 p~ ai '-> u 1 {v)=p' 2 ^ + (i- P - 1 ) 2 p~ 2a1 ' ( 8 ) 

7'=7+l 7'=7+l 

The integral in ([7]) can be easily evaluated. Keeping only the leading (quadratic) terms of expansion in a, we get 

Q{\ 7) = ^T+^ exp \- l -^f) (9) 

The approximation obtained for the expansion is valid for (^(S,^) 2 ^ <C 1, i.e. for a 2 <C 1. Let us note that 

it would be interesting to continue the expansions in ([3]) and in ([7]) keeping explicitly the terms of order of er 4 . The 
computation of the spectral density in this case involves the Airy functions. The discussion of this question will be 
published in the extended version of our work. 

Since the distribution function P(^ 7 n ' > ) docs not depend on n, the expressions (UJ-© hold for any of 1 < n < p r ~ 7 
eigenvalues on the hierarchical level 7. Introducing the number of eigenvalues on the level 7 (i.e. the degeneracy), 
.9(7) = P r_7 j we can rewrite the spectral density, p(X), now as follows 



1 r 



7= 



We have computed numerically p(\) for the Gaussian distribution -P(£ 7 "^), setting explicitly = p ' a+1 ' 7 e ^ . 
Simultaneously we have considered the linearized case of matrix elements, q 7 "' ) ~ p - ( Q+1 ) 7 (l — ^ 7 ™' ) ). We have 
compared in Figf2^ the explicit and linearized cases for T = 8 and a = 0.1. The same results for a = 0.5 are shown 
in FigfSb- We see that for a < 0.2 the linearization of the matrix elements suggested in ([3|) still makes sense for 
Gaussian distribution of energies. The analytic results computed using ([9|)- (|10j) coincide with the linearized case and 
are indistinguishable in the figures. 




FIG. 2: (Color online) Spectral density p(A) for V = 8: (a) a = 0.1: solid line - explicit case and dotted line - linearized case 
for a — 0.2; dashed line - explicit case and dotted-dashed line - linearized case for a — 0.5; (b) a = 0.5: notations the same 
as in Fig[2h. 



The analysis of pTT)) for arbitrary a > in the limit r — > 00 allows one to extract the asymptotic behavior of the 
tail of the spectral density p(X) as A — > —00: 



P(A) 



A— > — 00 




for a > 1 
for < a < 1 



(11) 
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The information about the relaxation characteristics of any system (including the hierarchical one) is encoded into 

the so-called "survival probability", W(t), having a sense of the probability to find a dynamical system in the initial 

state by the time t. For a nonrandom kinetic Parisi matrix of T hierarchy levels, the survival probability is simply 
r 

W(t, r) = (p — 1) P^ 1 eXjt +P _I \ where the eigenvalues, A 7 are defined in fT} (recall that A 7 < for 7 = 1, T). 
. 7=1 

For LCP transition matrices randomly distributed over the ensemble, the survival probability is a random quantity. 
To make the survival probability independent on the specific initial point on the tree, we consider W(t, T) averaged 
over all positions of initial states: 

r 

W(t, T) = ( P -1)J2p^ + = (P - 1) P^ 2 ^' 1> r ) + P~ T ( 12 ) 

{7,™} 7=1 

where 

Z(t, 7j r)=^e^- t (13) 

n=l 

The sum in (p~3|) runs over all p r ~ 7 tree vertices of the level 7 and the eigenvalues, A 7iH , are defined in ([4]). Note that 
despite W(t,T) does not depend on particular initial states on the tree, it is still a random quantity since Z(£, 7,r) 
and, hence, W(t,T) depends (via Eq. ([4j)) on the ensemble of random barriers, {£ 7 "' > }. 

The direct computation of (W(t, T)) is a difficult problem, though some important information about the proba- 
bilistic behavior of (W(t, T)) is available by analyzing the distribution functions of the sub-sums Z(t, 7, T). Note that 
the eigenvalue, A 7j „, of the kinetic LCP matrix T defines the escape rate from the "basin" (i.e. sub-tree) with the 
root located in the point ( 7 ,n) - see FigfTp. Hence, 7, T) is the partition function of all escape rates at time t 
from the hierarchy level 7. 

Further consideration is based on the observation that Z(t,j,T) resembles the partition function of the directed 
polymer on the disordered tree (DPDT) analyzed in Substituting the expansion (J3j) into (fl"3|) . and keeping 

only the linear terms in £ 7 ,l) (|£ 7 n) | < 1), we get Z(t,^f,T) = Z(t, 7, r)e"^ (r) *, where v 7 (T) is defined in © and is 
independent on any particular path on the tree, while Z(t,"f,T) has the following expression: 

Z(t, 1 ,T)= ux p| tL-^ + il-p- 1 ) P~ ai 'tf ] ) \ (14) 

all paths [ V 7 '= 7 +l / J 

The first summation in (|14|) is taken over all p r_7 vertices of 7-level on the p-adic tree, and the sum in the exponent 
runs along the path on the tree from the vertex (7, n) towards the root point (exactly as in ((3])). 

The function Z(t,j,T) satisfies the stochastic recursion (compare to [Tij ]) 

Z(T) = cxp [t(l - p-^p-^r] E Z 3 (T 1) (15) 

i=i 

Introduce as in [3] the averaged characteristic function, G m {x) by 

G m (x) = (cxp{-Z(m) cxp [til-p-^p-^x]}) (16) 
where 7 < m < T. Using the factorization of G m {x) on the Cayley tree, we come to the recursion for G m {x): 

G m -i(x) = J d$P(0 [G m (p- a (x + 0)] P (17) 

where P(^) is the Gaussian distribution ([6]) for £. To equip p7|) by the initial condition, put in (fi"4| 7 = T and get 
Z(t,-y = r,T) = e tp ~" « r . Thus, we fix the initial condition at the root of the tree by: 

G r (x) = <fflq){-fflq)[*(l-p- 1 )p- 0,r x] exp [tp- aT ^] }) p(fr) (18) 
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The linear approximation in (|14[) implies the narrow distribution P(£) with a 1 <C 1. In this case the boundary 
condition (fT5|) reads 

G r (x) = exp {- exp [t(l - p- 1 )p-° tT x\ } (19) 

Eqs. (fT7|) . Q19[) set the problem. It is more convenient to revert the direction along a tree by defining n = T — m 
(0 < n < T — 7). The case a > is a "contracting map" and will be analyzed elsewhere, while here we pay 

attention to the case a = (i.e. for (<Zt* / = P~ J )i which is formally identical to the Dcrrida-Spohn model of 

DPDT It is known from [li | that the solution of (17]), dT9j) in the continuum approximation is related to the 

Kolmogorov-Petrovsky-Piscounov (KPP) equation [l5[ and has a travelling wave. 

In terms of the work [l4j|, for long distances on the Cayley tree, the partition function, G n (x), is a travelling wave of 
the form G n {x) = w(x— fn), where the speed / is fixed by the initial condition, i.e. by the "inverse temperature", /3, 
where f3 = t(l— p^ 1 ) (for a = 0). The speed, /, is the free energy of the system per unit length in the long length limit 
f{t) = i In (pj d£,P{£)e K ). For t < t CI the travelling wave propagates with the speed f(t), while for t > t CT the speed 

/ is freezed at the critical value f{t cr ). The value i cr is determined by the solution of the equation ^/M| t _ t = 0. 
The overlap between two trajectories starting from the common root point can be: i) (with probability 1) for t < t CI , 
and ii) either (with probability ir(t) = const/i) or 1 (with probability 1 — 7r(i)) for t > t cr . 

Despite the formal correspondence of our problem for a = with DPDT is complete, the interpretation of the 
behavior observed in [l4| in terms of our model is not straightforward and deserves some discussion. The principal 
difference between the partition function (|14|) and the one of DPDT consists in the following. Our kinetic problem is 
defined on an p-adic tree, i.e. only on the boundary of an ultramctric p-adic the Cayley tree and all kinetic properties 
of our model deal with probabilities to cover some distance along a tree 's boundary for given time interval. To the 
contrary, the DPDT model is defined in the "bulk" of the Cayley tree, i.e. on the whole set of tree's generations, m, 
and the travelling wave propagates along m, from the root point, m = 0, "downwards". 

The overlap of the trajectories computed in [l4| allows one to connect the "bulk" and the "boundary" behaviors. 
Since on small time scales t < t cr (i.e. for "high temperatures") the overlap of any two trajectories (starting from the 
root of the tree) is zero, the boundary of the Cayley tree is uniform since there is no clustering. Hence, the escape 
probability from any basin at the boundary of the Cayley tree, is independent on the specific point on this boundary 
and the kinetics happens as for nonrandom Parisi matrix. To the contrary, on large time scales t > t CI (i.e. for "low 
temperatures") the overlap of the trajectories signals the "probabilistic non- uniformity" of the tree boundary: with 
the probability n(t) the boundary is still uniform, while with the probability 1 — Tr(t) the boundary has a "lacunary" 
structure such that the escape from some boundary basins is locked and the kinetics is confined in these basin for 
infinitely long time. 

The linear approximation ([3]) applied to (| 14[) deserves some special care. Actually, due to ([3]) the unphysical positive 
tail in the spectral density p(X) could appear. Since the survival probability W(t, T) depends on all eigenvalues (see 
Eq . (fT!?)) ) . for sufficiently large t these unphysical tail could have a non-negligible effect in the observed behavior. In 
order to eliminate such an influence, let us note that as a 2 decreases, the linearized behavior becomes better (see 
Fig|2]). So, we fix a and define the maximal time, i m ax(c), until which our predictions are valid and the interference 
of unphysical tail is negligible. The estimation of t m ax(c) is very straightforward. Namely, it follows from (fT2|) that 

the major contribution to W(t,T) at large t comes from the relaxation modes, for which t ~ 1. At the same 
time one sees from ([§]) that the unphysical tail corresponding to A > is exponentially small for w 7 3> a^Ju^. For 
practical purposes we may use the estimate v-y > 2a^Ju^j. Using ([8]) we get a < 0.5. Thus, for a < 0.5 the influence 
of unphysical tail is exponentially suppressed for any t. 

In this letter we have analyzed the simplest statistical properties of randomized ParisHype kinetic matrices. We 
also believe that the observed relation between dynamics on ultrametric landscapes described by LCP matrices and 
directed polymers on disordered trees is very promising for establishing a connection between the elements of p-adic 
analysis and disordered models on tree-like structures. Our work should be considered a preliminary step in this 
direction. The straightforward generalization of the obtained results concern an explicit account for the terms of 
order of Ifr™'*! 2 in ((3J) for the spectral density and for the survival probability. 

We are grateful to K. Bashcvoy and O. Vasilycv for illuminating numerical simulations, to O. Bohigas and Y. 
Fyodorov for helpful discussions, and to M. Mczard for paying our attention to the relation between LCP systems 
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